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Abstract 

We examine charged slowly rotating perfect fluids in the presence of a 
cosmological constant. The asymptotic form of the vacuum solutions to 
the linearised Einstein-Maxwell field equations is found and the possibility 
of matching this vacuum to the slow rotating Garcia metric is considered. 
We show that, contrary to the case of zero cosmological constant, this 
Garcia metric can be matched to an asymptotically de Sitter vacuum in 
the slow rotation limit. We conclude the Garcia metric may potentially be 
suitable for describing a charged isolated rotating body in a cosmological 
background. 


1 Introduction 

It is remarkable that one hundred years after the general theory of relativity 
was first formulated that there are still no known exact analytical solutions to 
Einstein’s equations describing an isolated rotating body. There are of course 
highly accurate numerical solutions of rotating stars mm, but mathematically 
it is of interest to obtain analytic solutions. In particular, and of main inter¬ 
est in this paper, there are no exact solutions describing an interior rotating 
charged perfect fluid which can be matched to an asymptotically flat vacuum 
exterior. One reason for this is the sparsity of known analytic charged rotat¬ 
ing fluid interiors. One such interior is the Garcia solution [3], also known as 
the Wahlquist-Newman metric [4], which is the charged generalisation of the 
rotating Wahlquist solution The Wahlquist solution is in turn the rotating 
generalisation of the static Whittaker metric. These metrics all have unphysi¬ 
cal equations of state, in the Wahlquist and Whittaker case the simple relation 
^ + ip = const holds, where p is the energy density and p is the pressure. 
Nonetheless due to the lack of explicit solutions the possibility of matching 
these solutions to external vacuum domains should be investigated. 
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The slow rotation formalism developed by Hartle in 1967 [B] has been impor¬ 
tant in deriving results about the possibility of matching these interior solutions 
to exterior vacuum regions. If the matching is not possible to a particular ex¬ 
terior in the slow rotation limit, then it will not be possible for a more rapid 
rotation. This argument has been used to show that the Wahlquist metric can¬ 
not be matched to an asymptotically flat vacuum [7], and one proves this by 
expanding the Wahlquist metric to second order in the fluid’s angular velocity. 
If one drops the requirement that the vacuum is asymptotically flat, the match¬ 
ing becomes possible [5] , however this means the Wahlquist metric cannot serve 
as a model for an isolated body, as a far away quadrupole mass distribution is 
required to keep the body in equilibrium. 

The matching conditions become more restrictive once charge is included. 
In the absence of charge, if one linearises the Einstein equations in the angular 
velocity of the metric, one can show that the general slowly rotating exterior 
vacuum solution of a perfect fluid sphere coincides with the Kerr solution, and 
one can match any slow rotating perfect fluid interior to this. However in 
the presence of charge, the general exterior vacuum solution to the linearised 
Einstein Maxwell equations does not in general coincide with the Kerr-Newman 
solution, and in fact is in general not asymptotically flat [S] . Thus if one requires 
asymptotic flatness, as necessary for describing an isolated rotating body in a flat 
background, then the matching conditions become overdetermined in general. 
In the particular case of the Garcfa metric, it was shown in [S] that to first 
order in angular velocity this metric cannot be matched to an asymptotically 
flat exterior metric. It is found the fluid needs to be embedded in an external 
magnetic field, parallel to the axes of rotation, and thus the authors conclude 
that the Garcfa metric cannot be suitable for describing an isolated charged 
body. 

It is now consensus that the universe is accelerating in expansion, which 
can be well explained by the presence of a small non-zero positive cosmological 
constant [iniin]. Thus the universe is not asymptotically flat, rather asymptot¬ 
ically it behaves like de Sitter space. This observation has resurrected interest 
in studying such spacetimes, which are either asymptotically de Sitter or anti- 
de Sitter. Of course, the small observed value of the cosmological constant 
means that its effects on astrophysical objects are negligible. However the pres¬ 
ence of even a tiny cosmological constant completely changes the asymptotic 
structure of spacetime [mill]. Moreover a given perfect fluid solution to Ein¬ 
stein’s equation in the absence of a cosmological constant is also a solution to 
Einstein’s equation with cosmological constant, achieved by making the substi¬ 
tution p —>• p —A/{8nG) for the pressure p, and p p +A/{SnG) for the energy 
density. This of course changes the equation of state of the fluid, but for the 
Wahlquist metric, the equation of state retains the same form, p + ip = const. 
Therefore the possibility of matching rotating fluids to asymptotically (anti)-de 
Sitter spaces should be investigated. Recently this approach was considered 
in |I2) . where it was shown that to second order in the angular velocity the 
Wahlquist metric can be matched to an asymptotically (anti)-de Sitter exter¬ 
nal vacuum. The purpose of this paper is to investigate this possibility for the 


2 


Garcia metric in the slow rotation limit. 

This paper is organised as follows. We begin in Section 2 by reviewing the 
Garcia solution and its form in the slow rotation limit. In Section 3 we examine 
the general slow rotating electro-lambda-vacuum exterior metric, which we show 
is always asymptotically (anti)-de Sitter, with the electromagnetic field tensor 
likewise asymptotically decaying. Finally in Section 4 we describe the matching 
procedure and we show that up to first order in the angular velocity one can in 
fact match the Garcia metric to an asymptotically (anti)-de Sitter vacuum. Thus 
we conclude by noting that the Garcia metric may be suitable for describing an 
isolated rotating charged body in (anti)-de Sitter spaces. 

2 The Garcia solution 

In this section we briefly review the Garcia solution and dicuss its form in the 
slow rotation limit. The Garcia solution, first discovered in 1991 [3], is the 
charged generalisation of the Wahlquist metric [S] and is given by 


^(dr -I- SNda)^ - + SMdaf' -I- -|- 


where 


A = M-N, 



In the above a, &, e, g, k, m, n, a, /3, 5 and are eleven constants. Up to diffeo- 
morphism invariance only eight of these are independent parameters [9] . 

This metric is a solution to the Einstein-Maxwell equations with a cosmo¬ 
logical constant 


Gap + ^QaP — ^'^{Tap + ^ap)^ 

VpF^^ = J“, V[aFp^] = 0, 


( 2 . 1 ) 

( 2 . 2 ) 


where is the electromagnetic field tensor which can be written in terms of 
the four potential as 


F^1/ — ^ jyA^. 


(2.3) 


is the stress energy tensor of the electromagnetic field 



(2.4) 
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(2.5) 


and is the stress energy tensor of a standard perfect fluid with 

'^af3 = (M + P)UaU/3 + pQaP, 


where is the four velocity. 

The pressure p and density p of this particular solution are given in terms 
of the metric components by 



( 2 . 6 ) 


(2.7) 


where the quantity S has been defined as 



( 2 . 8 ) 


The four potential A is given by 



(2.9) 


( 2 . 10 ) 


Now we examine its form in the slow rotation limit. One can write the 
general metric of a slowly rotating perfect fluid sphere in the following form [5] 


+ C^[d6^ + sin^ 9{d(j) — wdt)^] (2-11) 


where an expansion in the angular velocity VL of the fluid is made. In this 
paper we consider only the first order perturbations in 17. This means that the 
metric functions A, B and C are simply the static non-rotating solutions, to first 
order in VL they remain unmodified. Rotation is added through the function lo 
which is assumed to be of 0{Vl) (and thus we follow others and slightly abuse 
notation by neglecting the component of the above metric). In coordinate 
transformations of the Garcfa metric were derived which allowed it to be put 
into the above slowly rotating form. It was found that 



where the functions hi and w are given by 



(2.13) 



sin z 


(2.14) 
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and ro is of order and is taken as the slow rotation parameter. In the static 
limit To —>■ 0 this metric reduces to the charged Whittaker metric. The density 
and pressure are now given by 

8ttp = + ^ - /3^) + 4/37^At + A 

Snfi = 72 ( 3/11 - 4 + /32) - 4 / 372^4 - A. (2.15) 

We must also expand the four potential A up to first order in rg, and it is found 
that 

A = —j3z cot zdt + rg [/3 sin2 0(1 — 2 : cot z) — (i — g cos 9]dip. 

In this expression g is the magnetic monopole charge parameter. The metric 
is completely independent of this parameter, and it only contributes to the 
Maxwell equations. 

The set of constants in the original metric have now been transformed to a 
new set of five constants /3, k, 7 , g and rg. The number of constants has been re¬ 
duced from the original eight independent parameters (modulo diffeomorphisms) 
by fixing a coordinate system and ensuring that the metric is completely regular 
at the center. 


3 Electrovacuum exterior 


Now in this section we will find the corresponding slow rotating form of the 
vacuum exterior. We make a slow rotation approximation around Reissner- 
Nordstrom-de Sitter space by expanding the metric up to first order in the 
angular velocity fl. We write the metric in the same form as the slow rotating 
Garda metric ()2.11|1 . Now from the Einstein equations, to first order in the 
angular velocity, the diagonal components of the metric are unmodified and the 
only modification is to the gtip component, see [ 6 ]. Hence we may write the 
electrovacuum exterior as 


, 2 , 2M e2 Ar2 „ 

ds = — (1 - 1 —^ - —)dt 

r r 0 

+ r2[d02 + sin2 9{d(p — wd/)2], 


2M e2 Ar2 ^ 

(1-+ ) idr2 

r 6 


( 3 . 1 ) 


where M is the mass, e is the charge and A is the cosmological constant, which 
we consider either to be positive or negative. The function w is first order in the 
angular velocity, and again we ignore the ^2 component of the above metric. In 
the Garda interior solution, the function w is dependent on r only, and we will 
make the additional assumption that this to be true in the exterior metric also. 

We can choose a gauge such that the four potential of the electromagnetic 
field Aa = Aa (r, 9 ), to zero-th order in angular velocity, is 

At = —, Ar = 0 , Ag = 0 , A^p = 0 . ( 3 - 2 ) 

From the Einstein equations we will see that the only component of the four 
potential that is perturbed to first order in H is A^. 
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3.1 Slow rotating Kerr-Newman-de Sitter metric 

Let us first examine the slow rotating form of the Kerr-Newman-de Sitter metric. 
The full metric in Boyer-Lindquist coordinates is given by 


ds^ = — asin^ -|- — (r^ -|- a^)d(j)Y 


where 


-I Pi 


-I Pi 


+ (3-3) 

A/ fc)/ 


pi = \/r"^ + a? cos^ 9, 


(3.4) 

0; = 1 + yy cos^ 9, 


(3.5) 

r"2 

u = --pV + a ) + r 

- 2Mr + a^ + e^. 

(3.6) 


And the electromagnetic potential of this metric is given by 


A = - 


er 

pI 


sin^ 0 ^ 

dt — a——dip 


(3.7) 


To write this in the slow rotating form (12.111) , we simply expand this metric 
up to first order in the parameter a, which we take to be the rotation parameter 
and is of order 12. It is readily seen that the functions and w are given by 


.KN 



(3.8) 


UJ 


KN 


2aM ae^ aA 

J.3 3 


(3.9) 


We are free to perform a rigid rotation of our coordinate system, p ^ p + cit, 
where ci is a suitably chosen constant to remove the contribution a A/3 from 
this expression. This is done so that we are working in an asymptotically non¬ 
rotating frame. 


3.2 General vacuum 


Let us now return to working with a general exterior. For a general A-electrovacuum, 
the {p, t) components of Einstein’s equation gives us the following differential 
equation 


. -iduj 4e dA^, 
dr^ dr sm 


9 dr 


= 0 . 


(3.10) 
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We can solve this equation to find the general solution of A^p for a non-zero 
charge 

+ f{9), (3.11) 

where / is an arbitrary function of 9. 

We now insert this solution back into the Maxwell equation = 0 and 

we find the following differential equation relating w to / 


— {2Mr — H--|- {-Mr — 2r^ — -e^ 

4 ^ 3 ' dr^ ^2 2 


5A 4 ^d?u) 

2 r >’■ ^ 


+ (4Ar + y-r = ~^a 

2 6 dr sin 9 




df{9) 

d9 


)■ 


(3.12) 


This is a separable equation, with the left-hand side independent of 0, and hence 
we may write 


1 dV( 0 ) 


— cot 9 


df{9) 

d9 


) 


K 


(3.13) 


where is a constant. We can write the left-hand side of this equation simply 
as 


1 ( d'^f{9) ^^.p dfjd) . ^ d^f {cos 9) 

sin^0 d9^ d9 d{cos9)'^ ’ 

so that this equation (13.131) has the simple general solution 

lATcos^e Ci ^ 

f\d) = o-^-^-cos6». 

2 e e e 


(3.14) 


(3.15) 


Now examining the radial part of the equation (13.121) gives us a third order 
inhomogeneous differential equation for ui 


{2Mr — — r^ + + {i^AIr — 2r^ — —e~ 


id^U! 

dr^ 


- r^)r^ - 

2 3 ’ dr^ 


/,,, 5 -idu} 

+ {AM — -r + 5—r^)r^— = K. 
2 6 dr 


(3.16) 


The presence of A in this equation makes finding an analytic solution difficult. 
However we can find a particular solution, and noting that a constant is also 
clearly a solution allows us to write the following ansatz for w 


uj = Cq + K 


\3Mr^ 



+ CiUJi{r) + C2U}2{r), 


(3.17) 


where Co, Ci and C 2 are all first order constants in H and wi and UJ 2 are 
two linearly independent solutions of (13.161) . We will analyse the form of these 
solutions shortly. Now since we know that the slow rotating Kerr-Newman-de 
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Sitter is a particular solution to the above equation, we see that we can set 
K = —3aM, so that taking Ci = 0, C 2 = 0 we have simply the slow rotating 
Kerr-Newmann-de Sitter metric. And again, we are always free to set Cq to 
zero by adjusting to an asymptotically non rotating frame via a rigid rotation 
of our coordinate system. 

Substituting this solution back into the potential we obtain 

= — sin^ 9 — gin^ 9{CiUj[ + C 2 OJ 2 ) + — + — cos0. (3.18) 

r 2e 4e e e 

The term Cg simply gives rise to a magnetic monopole. Now the form of the 
solutions uji and uj 2 of (13.161) are completely different depending on whether 
A = 0. We investigate the two distinct cases below. 


Zero cosmological constant A = 0 


The solution for A = 0 was derived in [9], and we review its properties here. 
Assuming ^ M^, it was found 

e^(r^ — 2Mr + e^)[r + M + (e^ — r^)L(r)] 2e^(e^ — 2Mr) 

“^1 =- TTTrys —- 1 " ( 3 . 19 ) 


M 2 (e 2 - M2)2r4 

2 e^(e^ — 2Mr) 

^2 = - + 

where L(r) is implicitly defined through 

dL{r) 1 


3M3r4(e2 - M 2 ) 


dr 2Mr — — e2 ’ 


(3.20) 


(3.21) 


and can be written in terms of log or arctan depending on whether the discrim¬ 
inant of the metric component gtt, which is given by 4(M^ — e^), is positive or 
negative respectively. Asymptotically expanding these functions in powers of 
1 /r we find the following far field behaviour 


2 e^ 4e^ 

^ 3M3r4 5M2r3 

2 2e^ 

W2 —-15—h . -|- . . . 

r r3 Mr^ 


(3.22) 

(3.23) 


Now we can find the asymptotic behaviour of the metric and electromagnetic 
held tensor by using the following tetrad 


eo = (1 - 


2M 


+ f!)-i/ 2 ^ 

r2 dt 


dr 

1 d 
“ rm 

Id . 2M 

63 = —r-x ^-wr sm 0(1- 

r sin u Oip r 




dt' 


(3.24) 














Expressed in this tetrad, it can be shown the electromagnetic field does not 
vanish asymptotically unless the constant C 2 vanishes, since 

lim Fi 3 = -^sin 6 >, (3.25) 

r—^oo c 

lim F 23 = — —cos 9, (3.26) 

r—>-oo g 

and similarly by expanding the Weyl curvature tensor in this tetrad, it is also 
found that the metric is not asymptotically flat unless C 2 = 0. Thus the general 
vacuum exterior rotating charged metric cannot describe an isolated object. It 
was shown in that if one attempts to match the Garcia metric to this exterior 
vacuum, the matching conditions dictate that this constant C 2 7 ^ 0 , and so the 
exterior vacuum will not be asymptotically fiat. Therefore the Garcia metric 
cannot describe an isolated rotating body in an asymptotically fiat spacetime. 


Non-zero comological constant A 7 ^ 0 

When A ^ 0 we are unable to find the exact analytic solutions for wi and 002 
from equation (13.161) . However we are able to find the form of the solutions 
asymptotically by making an expansion in powers of 1 /r 


1 6M 

- 7X^7 +- 

_ 1 9 3e2 9 

2 .^3 5Ar® 7Ar^ 7A^r’^ 


(3.27) 

(3.28) 


and we note the function uji now falls off like l/r'^ and UJ 2 now falls of like 1 /r^. 

To confirm these results, we can also integrate (13.161) numerically. We set 
AT = 0 in the equation so we ignore the particular solution corresponding to the 
Kerr-Newman-de Sitter contribution and find simply the form of the wi and W 2 
contribution. We also numerically integrate the solution for uj'(r), since we are 
only interested in how the solution decays for large r, and this will allow us to 
ignore the constant contribution which arises when numerically integrating w. 

A plot of uj'{r) against r is plotted on a logarithmic scale in Fig. ([TJ for 
a negative cosmological constant. A positive cosmological constant leads to 
numerical difficulties at the cosmological horizon, meaning different coordinates 
would need to be used, however numerically integrating the equation above the 
horizon gives very much the same result. On logarithmic graphs, the functions 
r" are straight lines with gradient n, so analysing the gradient of the line on the 
logarithmic scale allows us to measure the power at which r decays. We find for 
all three graphs n ~ —3.9, which supports our asymptotic result that uj' decays 
like \/t~^ for large r. 

Now we must investigate whether this observed asymptotic behaviour means 
that the metric is asymptotically (anti)-de Sitter, and we must also investigate 
the asymptotic behaviour of the electromagnetic field tensor. This time we use 
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Figure 1: Plot of u>'{r) against r on a logarithmic scale on both axes. The 
parameter values are A = —0.001 and e = 0.1 and M = 1. This means the 
black hole horizon is located at rbh ~ 1.99. The initial conditions are w(2) = 
{1,0.1,0.01}, u;'(2) = {-0.01,0.02,0.01}. 


the tetrad 


eo = |1 - 


2M 


( 2M 

1 d 

“ rm 

1 

63 = 


d 


Ar^ 

Ar^ 


- 1/2 


dt 


1/2 


dr 


■ n,. 2 M 

. - wr sin 0(1-h — 

r sm 0 oip r 


Ar^ 1 d 
3 ^ dt 


(3.29) 


In this tetrad, we obtain the following behaviour of the electromagnetic field 
asymptotically 


lim Fi 3 = 0, (3.30) 

r—^oo 

lim F 23 = 0, (3.31) 

r—¥oo 

and in fact these components decay like 1/r^, exactly the same as Kerr-Newman- 
de Sitter space. This is in stark contrast to the A = 0 case. 

It is also readily seen that this exterior metric is always asymptotically (anti)- 
de Sitter. We can see this by examining the asymptotic behaviour of the Weyl 
tensor. We use the standard definition for the gravitoelectric fields Ei and the 
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gravitomagnetic fields Hi 

El = i?0101 ~ ’ E 2 = i?0202 ~ j E 3 = i?0102, (3.32) 

Hi = i?0123, H 2 = —i?0213, H 3 , = i?0223, (3.33) 

with the remaining four components of the Weyl tensor vanishing identically. 
Using the asymptotic expansion for w, we find the following expansions for 
large r 

E2=0{r-\n^), E3=0{Sf). (3.34) 


For the gravitomagnetic fields Hi it is found 


Hi 

H 3 


3^2 cos d 5 

2r4 ^ ’ 




H 2 = 


3 C 2 cos 6 


V^Czsinff 2 v^Ci sin 9 


2r^ 




2r4 


+ oir-^,n^), 


(3.35) 


which decays appropriately quickly to be asymptotically (anti)-de Sitter. The re¬ 
maining requirements of being asymptotically (anti)-de Sitter are easily checked 
and can be done in a similar way to the examples in [laiiaiii]. In the case 
of de Sitter space it is of course important to be careful with coordinates above 
the cosmological horizon, and this can be dealt with in a rigorous way, again 
see M- It is also possible to observe that the extra contributions to uj decay 
as quickly as Kerr-Newman-de Sitter space, which is asymptotically de Sitter. 

Thus the general vacuum exterior is asymptotically (anti)-de Sitter, and 
hence the possibility of matching the slow rotating Garcia metric to this exterior 
vacuum should be investigated, considering we now have one extra constant 
available to perform the matching with. We investigate this in the next section. 


4 Matching 

In this section we describe the matching procedure and prove the following 
result: 

The linearised Garcia metric can be matched to an asymptotically (anti)-de 
Sitter vacuum exterior. 

Now we have already written both the interior and exterior metrics in the 
form 

ds^ = -A^dt^ E B^dr^ + r^[de^ + sin^ 9(d(l> - ujdtf]. (4.1) 

In what follows we will denote the vacuum region and the interior region by 
superscripts (U) and (/) respectively. We first apply the following coordinate 
transformations to the interior fluid region 

t = Get' tp = If' + VLt (4.2) 
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in order to perform the matching. 

The matching conditions up to first order in the angular velocity of the fluid 
were derived in [^. Up to this order the matching surface, or the zero pressure 
surface, remains a sphere, which we take to be at radius r = ri. Deviations 
from the spherical symmetry of this surface don’t appear until one considers 
second order corrections to the metric. In order to have a global model in which 
surface layers of matter are absent, the following Darmois-Israel m conditions 
on the metric are required 


(V) (I) 

9ab =9ab 

K^v) _ ^(/) 
^ab — ^ab 


(4.3) 

(4.4) 


where Kab is the extrinsic curvature of the metric. These conditions imply the 
metric functions A and lo are continuous and B is continuous at the zero 
pressure surface r = ri. To ensure absence of surface charges and currents the 
components Frt, Fr^p and Fg^ of the electromagnetic field tensor are required 
to be C° on this surface. 


4.1 Zero-th order matching 

The zero-th order matching conditions allow us to find the constants Cq, M, e 
and the radius ri of the zero pressure surface in terms of the parameters /c, 7 
and /3 of the internal static fluid. Matching of the ggg component gives 


sin zi 

ri = - 

7 

Using continuity of Kgg and gtt we obtain 


Cq = cos zi 


and continuity of At^r gives 


e = 


— (zi — sin zi cos zi ). 
7 


(4.5) 


(4.6) 


(4.7) 


The equation of the zero pressure surface, which can equivalently be derived 
from the junction condition for gtt and Ktt, lets us solve for the constant 


1 A ^tanzi ^ 2 


7:7 = (1-77)- 


-|- 2/3^ (2 -I- zi cot 2zi). 


7 ^ Zi 

And finally the mass of the fluid is derived from the continuity of gtt 
2 7 ^ 3 

1 

H- (z? zi sin 2zi cos 2zi -\— sin^ 2zi ). 

27 sin zi 2 


(4.8) 


(4.9) 
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4.2 First order matching 

Now if we consider the first order matching conditions, we are able to find the 
values of the first order constants a, il, Ci, C 2 and in terms of the parameters 
of the interior fluid. Continuity of the gt^, component of the metric gives the 
angular velocity of the fluid 


n = —(/3^Zj(cot^ zi — 1 ) — /3^ + K ^{1 

sin zi 

, 2aMj^ 

H-(Ciw^ri) + C 2 a; 2 ri) + ^-3 - 

cos zi sin zi 


— zi cot Zi)) 

2 4 

ae 7 

“T^ ’ 

sin zi 


(4.10) 


while continuity of Ktip allows us to find the Kerr parameter a 


/-a 2 /■r>o2 2 2 J.3 

a = -——-r(7 >0 CSC zi sec zi (2/3 k z^ coC zi 

K^7^(4e^7 — 6M sinzi) 

+ cotzi(—2/3 ^k^(zi + 1) + 20^k^z\ csc^ zi + 3) — 2 zi(/3^k^ + 1) cot^ Zi 
+ Zi(2/3^k^ - csc^ zi)) - C'ia;i(ri) - C'2a;2(ri)). (4-11) 


Equating magnetic monopole terms gives, as in the A = 0 case, we find simply 


Cs = -roeg. 


(4.12) 


Finally the continuity of F^pr and F^g allow us to solve for both Ci and C 2 


Cl 


— -fi (<^2('Ci) (Sae + e7^r?ro csc^ zi seczi 

Di'° 

+ riuj 2 (ri)(ae - /3riro + /3riroZi cot zi)) 


4/3riro(l - zi cotzi)) 

(4.13) 


C 2 = (5«e + csc^ zi seczi — 4/3riro(l — Zi cotzi)) 

+ riw"(ri)(ae - l3riro + /3riroZi cot zi)) (4.14) 

where the quantity D is dehned as 

D = uj[{ri)uj2{ri) - w"(ri)w2(r'i). (4.15) 

We can then insert the solution for a (j4.1ip into (I4.13p and (I4.14L which will 
then give two linear simultaneous equations for Ci and C 2 which can then easily 
be solved, although the expressions are lengthy so we omit the details here. This 
means that all the matching conditions are now satisfied, and therefore our main 
result is shown to be true. 


5 Discussion 

We have shown that the Garcia metric can be matched to an exterior A-vacuum 
solution up to first order in the fluids angular velocity. This exterior solution is 
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asymptotically (anti)-de Sitter. This contrasts with the A = 0 case, where it is 
the case that the exterior of the slow rotating Garcia metric is not asymptot¬ 
ically flat in general. The addition of a cosmological constant essentially gives 
rise to an additional degree of freedom by allowing one extra constant in the 
exterior, meaning there are now enough parameters to match to all of the inte¬ 
rior parameters. This allows us to interpret the slow rotating Garcia solution as 
an isolated body in a cosmological background, an interpretation not possible 
in the absence of the cosmological term. 

This is a very similar result to the one recently obtained in [T^] for the 
Wahlquist metric, where to second order in angular velocity the matching to 
an asymptotically empty domain was found only to be possible once A was 
introduced. The addition of charge does not drastically change this result, 
however the requirement of needing a cosmological constant to perform the 
matching now appears at first order in the angular velocity rather than second 
order. In fact with this extra degree of freedom, it can easily be shown, using 
the same procedure outlined above, that any slow rotating charged perfect fluid 
interior can be matched to an asymptotically (anti)-de Sitter metric up to first 
order in angular velocity. 
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